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The experiments by Darbyshire and Mullin (J. Fluid Mech. 289, 83 (1995)) on the 
transition to turbulence in pipe flow show that there is no sharp border between initial 
conditions that trigger turbulence and those that do not. We here relate this behaviour 
to the possibility that the transition to turbulence is connected with the formation of 
a chaotic saddle in the phase space of the system. We quantify a sensitive dependence 
on initial conditions and find in a statistical analysis that in the transition region the 
distribution of turbulent lifetimes follows an exponential law. The characteristic mean 
lifetime of the distribution increases rapidly with Reynolds number and becomes inacces- 
sibly large for Reynolds numbers exceeding about 2200. Suitable experiments to further 
probe this concept are proposed. 



1. Introduction 

The transition to turbulence in pipe flow has been the subject of m any investigations 
since the first documentation of its phenomenology in l(Re vnoldsll 1 8 8 3l) . Reynolds already 
noticed that there was no sharp transition, that the transition could be delayed to very 
high Reynolds number when perturbations to the laminar profile were carefully avoided 
and that sometimes turbulence appeared without clear evidence for a perturbation that 
may have triggered it. His fi ndings have been confirmed and expanded on in many 
subsequent exper i ments, e.g. jWvgnanski fc Champagne! Il973t IWvgnanski et a/.lll975r 
iRubin et a/.lll980HEggels et a/.lll994HMa et aJl999HDarbvshire fc Mullinll 19951) On th e 
theoretical side, no linear ins tability could be found (see ISchmkTlTrfem'ihilrsorl l|l994f) : 
iMeseguer fe Trefethenl l|2003|) and references therein) but analysis of the consequences of 
the non-normality of the linearized problem has led to the identification of efficient ampli- 
fication me_chamsms_thjit caii^ perturbations of the laminar 
profile llBoberg fc Bros"alll988l ISchmid fc Henningsonlll99i iTrefethen et n/J 1200(1 Il993t 
iGrossmamfeOOftlHof et flil2003l). On the othe r hand the simulations in l|Brosall 989) and 
the experiments in ()Darbvshire fc Mullinll 995]) suggest that even if the perturbations are 
large enough to trigger turbulence, the flow can re-laminarize wit hout any previous indi- 
cation. Similar behaviour has been found in plane Couette flow llSchmiegel fc Eckhardtl 
1199711200(1 iBottin et o,/Jll99RHFaisst fc EckhardtlbflOOl: lEckhardt et «iJl2002l ~ 

Decay of the turbulent flow must follow from the dynamics of the fully developed 3-d 
turbulent state and cannot be explained by linearization around the laminar flow. The 
dynamical system concept compatible with such a behaviour is that of a chaotic saddle, 
a transient object with chaotic dynamics (jTeJ Il99ll IOtdfl993|) . The simplest example 
of a chaotic saddle arises for a particle in a box with a tiny hole: the dynamics inside 
the box can be chaotic, as measured by positive Lyapunov exponents, but it is transient 
and ends when the particle leaves through the hole. The escape through the hole is a 
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global event, and its rate is not related to the Lyapunov exponent. The analogy then 
is that the turbulent state is motion in the box, and escape from the box is related to 
relaminarization. The characteristic signatures of a chaotic saddle are: (a) a sensitive 
dependence of lifetimes on initial conditions, (b) an exponential distribution of lifetimes 
for initial conditions around the chaotic saddle, (c) a positive Lyapunov exponent in the 
chaotic phase, and (d) an independence of the variations of Lyapunov exponents and 
the escape rates with parameters. We will show here that the transition to turbulence in 
pipe flow shows all these characteristics and is thus compatible with the formation of a 
chaotic saddle. 

Our analysis is based on numerical simulations of the time evolution of perturbations 
in circular pipe flow with periodic boundary conditions in the downstream direction. 
The axial periodicity length will be too short to simulate localized turbulent struc - 
tures, such as puffs and slugs i Wvgnanski fc Champagnejll97.lt IWvgnanski et oiJll975h . 



but long enough to c apture the local turbulent structures. Following the simulations 
of feggels et a/Jll994) we take a length 10i?, with R the pipe radius. The turbulence 
will then typically fill the entire volume and we do not have to consider the advection of 
the turbulent state by the mean profile: lifetimes can be defined locally. The Reynolds 
number Re is based on the mean streamwise flow velocity U and the pipe diameter 2R, 

Re = 2RU/v. (1.1) 

and the units of time are 2R/U. As in the experiments l|Darbvshire fc Mullinlll99 *3) we 
keep the volume flux constant in time: this simplifies the analysis and prevents decay 
due to a reduction of flux that could occur in pressure driven situations when the flow 
becomes turbulent. 



2. Numerical considerations 

We use cylindrical coordinates and employ a pseudo-spectral method with Fourier 
modes in the periodic azim uthal and downstream direction and Legendre collocation ra- 
dially ijCanuto et aZjfl98^ . Various linear constraints on the velocity field are treated to- 
gether by the method of Lagrange multipliers: the rigid boundary condition, the solenoidal- 
ity, and the analyticity in the n eighbourhood of the coordinate singularity at the center 



line i Privmak fc Mivazakilll99sh . We verified our numerical scheme in various ways. For 



the linearized dynamics we reproduced the literature data on the eigenvalue spectrum of 
the li nearized Navier-Stokes operator with full spectral precision ijSchmid fc Henningsonl 
11994 . For the no n-normal linear dynamics and nonlinear 2D dynamics we reproduced 
Zikanov's results llZikanovl 1199^1 . At Re = 5000 a long, fully nonlinear 3D turbulent 
trajectory was analys ed and its statistical properties agreed with previous numerical and 
experimental results l|Eggels et q/Jll994t lOuadrio fc Sibillall2000h . 



The results presented here are obtained with azimuthal and axial resolution of |n/17| + 
|to/15| < 1, where n and m are the azimuthal and streamwise wavenumbers, respectively, 
and 50 Legendre polynomials radially. This resolution is not the maximal that could be 
integrated for a single run but reflects a compromise between the mutually exclusive 
requirements of maximal resolution, maximal integration time for a single run and a 
large number of sample runs for the statistical evaluations. It is justified by comparisons 
with lower and higher resolutions which show no significant differences for the range 
of Reynolds numbers Re < 2300 to which we restrict our computations. The maximal 
integration time after which we terminate a turbulent trajectory is 2000, sometimes 
3000 natural time units. This by far exceeds the values access ible in the longest currently 
available experimental set up l|Ma et «iJll999HHof et ali2QQ$f . For the statistical analysis 
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Figure 1. Traces of normalized pressure gradient (top) and energy content E^d (Eq. 12,21 
(bottom) for three slightly different initial conditions at Re — 2100. The initial energies are 
chosen high above the turbulent mean and the trajectories require about 200 time units to relax 
to the neighborhood of the turbulent mean. Pressure drop and energy content show correlated 
fluctuations. Note that within statistical fluctuations the averages of the turbulent dynamics 
are the same. This includes the dashed trajectory which happens to decay near t = 2000. The 
bottom frame shows the kinetic energy of the streamwise modulated part of the velocity field, 
Ezd, normalized by its turbulent mean as a measure of the velocity fluctuations 



and the calculations of Lyapunov exponents more than 1000 runs were needed, adding 
up to several years of CPU time on a single 2.2 GHz Pentium4/Xeon processor. 

The initial conditions for each run are the parabolic profile uhp to which a three- 
dimensional perturbation u p is added. The specific form of u p should not matter, as long 
as it triggers a transition to turbulence, since the dynamics of the turbulent state has 
positive Lyapunov exponents and leads to a quick elimination of any details and memory 
of the initial conditions. Furthermore, experiments with different k inds of wall-normal or 
azimu thal jets or other perturbations lead to very similar results l|Darbvshire fc Mullinl 
Il995lh so that we can safely assume that there is only one turbulent state to which all 
initial conditions are attracted. Therefore, we take as initial condition an uncorrelated 
random superposition of all available spectral modes. Changes in the initial conditions 
are then limited to variations in the amplitude of the perturbation, not in the form, i.e., 
throughout this paper we scan the behaviour along a one-dimensional subset in the space 
of all velocity fields, u(0) = uhp + a,u p . 

We define the amplitude A of an initial disturbance as its kinetic energy in units of 
the energy of the laminar Hagen-Poiseuille profile uhp, 

A = a 2 / / . (2.1) 
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The measures that we apply to the turbulent runs are the pressure gradient required 
to maintain a constant mean flux and the energy variable E^d, the kinetic energy in the 
streamwise modulated part of the velocity field, 
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E ™ = Y.Y. L |un, m (r)| 2 rdr. (2.2) 
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Figure 2. Turbulent lifetime vs. perturbation amplitude for different Reynolds numbers. For all 
Reynolds numbers a smooth region with short-lived states can be observed for small amplitudes. 
For larger Reynolds numbers this is followed by a ragged region of highly fluctuating lifetimes. 
The threshold amplitude that divides the two regions decreases with Reynolds number. The 
little rectangle at amplitudes of about 0.03 at Re — 2000 is further magnified in Fig. [3] 



The significance of is that if it becomes too small, then the flow is too clo se to 
an axially translation invariant flow field and will eventually decay (|Zikanovlll996|) . We 
therefore terminate integration if it drops below a threshold of 10 -4 . 



3. Turbulent lifetimes 

Typical trajectories are shown in Fig.^ Within about 150 time units they relax towards 
the turbulent state. The dashed trace belongs to an initial condition for which the lifetime 
of the turbulent state is about 2000: there is no indication of the decay until the energy 
E^d drops so low that turbulence cannot recover. 

Scanning the lifetimes of turbulent states as a function of initial amplitude for different 
Reynolds number we obtain the results shown in Fig. [2] For sufficiently small amplitudes 
all states decay and the lifetimes are short. Beginning with a Reynolds number of 1800 
sharp peaks with very long lifetimes appear. Beyond Re = 2000 several initial conditions 
reach lifetimes up to the integration cut-off of t max = 2000. Moreover, the variations in 
lifetimes between neighboring sampling points in amplitude increase. 

The transition between the short lifetimes for the initial conditions that decay quickly 
towards the laminar state and the longer ones for trajectories that show some turbulent 
behaviour is a very rapid one. Successive magnifications of a small region in amplitude 
for Re — 2000 are shown in Fig. [3 The increase is not monotonic, with modulations and 
structures superimposed. 

Fig. ^compares the energy traces for four initial conditions in two ranges of amplitude. 
Fig. 0Jl for the smaller mean amplitude shows that in an interval of width about 3% of 
the mean amplitude the lifetimes increase by about 40%. At a slightly higher amplitude 
Fig. shows that in an interval of relative width 3 • 10~ 4 the lifetimes increase by a 
factor of almost 4! 
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Figure 3. Successive magnifications of the lifetime variations with amplitude in the transition 
region for Re — 2000. The left panel is a magnification of the box indicated in Fig.|5J the right 
one a magnification of the interval indicated in the left panel. Note the increased ordinate scale 
in the right panel. 
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Figure 4. Energy traces for trajectories in small intervals of amplitude close to the transition 
region at Re = 2000, marked by open circles in the corresponding panel in Fig. [!|] The small 
differences in initial conditions are amplified and lead to noticeable differences in their evolution. 
For panel (b) the lifetimes are larger and the dynamics shows several large amplitude oscillations 
connected with the turbulent regeneration dynamics. 



The strong fluctuations in lifetimes are not limited to variations with amplitude alone, 
they also occur when changing the Reynolds number for fixed amplitude, as shown in 
Fig. [3 Such fractal behaviour in the variations of lifet imes under parameter variation 
has previously been observe d in plane Couette flow l|Schmiegel fc Eckhardtll997j) . Taylor- 
Couette flow for large radii ijFaisst fc Eckhardtl2000HEckhardt et all2002l). an d in models 
with a sufficient number of degrees of freedom l|Eckhardt fc Mersmannl ll999'). 



4. Lyapunov exponents 

The sensitive dependence on variations of initial conditions can also be quantified in 
terms of Lyapunov exponents. For a trajectory u(t) and deviations 5u(t) the largest 
Lyapunov exponent is defined as 

A= ton lim -log Ml (4.1) 

||«„(o)||-o t^oo T B ||<fa(0)|| V ; 

where | • 1 1 denotes the Euclidian norm. Instead of a single A from an infinite time interval 
we calculate ensembles of finite-time Lyapunov exponents, extracted from the integration 
of a small deviation from the reference trajectory. The separation was measured after 



1500 



500 



2%l 





I I I I 
(a) 


i j i | i 


TT 1 1 | 1 1 .T 








f; 










1? 


• 


















1 A; 
















i i i i 


- 

II 



Holger Faisst and Bruno Eckhardt 

2000 



2010 
Re 




500 



2067.0 2007.2 2007.4 2007.6 2007.8 2008.0 
Re 



Figure 5. Variations of turbulent lifetimes with Reynolds number for a fixed initial distur- 
bance amplitude A — 0.00327. The parameter range marked in (a) by the shaded interval near 
Re = 2007 is magnified in (b). 
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Figure 6. Lyapunov exponents in pipe flow, (a) Local Lyapunov exponent for a Reynolds 
number of Re = 2150. The fluctuations in the local Lyapunov exponent correlate strongly with 
the energy content E of the disturbance (in units of the energy of the laminar profile), (b) 
Largest Lyapunov exponent of the turbulent state with an error estimate based on the root 
mean square deviations of the ensemble of finite-time Lyapunov exponents. The line indicates 
a linear fit with slope dX/dRe ~ 4 x 10 -5 . 



time intervals of T — 200 units, and then rescaled to | |<5u| | = 10~ 9 . An initial time interval 
of 200 units was omitted because of the transient relaxation of the initial conditions onto 
the turbulent state. Similarly, the last 200 time units were omitted to avoid the decay to 
the laminar state. The infinite-time Lyapunov exponent was then approximated as the 
average of the finite-time Lyapunov exponents. 

The instantaneous changes in the sepa ration between neighbo uring trajectories can be 
used to define local Lyapunov exponents IjEckhardt fc Yaoll993|) . They correlate strongly 
with large energy fluctuations, see Fig.HJ.. When new large scale structures are generated 
the energy grows strongly and the Lyapunov exponent increases. Towards the end of a 
nonlinear regeneration cycle the energy goes down and the Lyapunov exponent decreases 
as well. Therefore, the fluctuations in the exponent are large and averages over at least 
10 time units are needed for the determination of reliable numbers. Obviously, this is 
very difficult to achieve for Reynolds numbers below Re — 2000, where only very few 
trajectories stay turbulent for sufficiently long times. 

The largest Lyapunov exponent is shown in Fig. 03. Its typical value is about 6.5 x 10~ 2 
at Reynolds numbers around 2000. For predictability this Lyapunov exponent implies 
that an uncertainty doubles after a time of (In2)/A « 10 units. Over a time interval 
of 200 units a perturbation grows 10 6 -fold! Note that 200 units is about the separation 
between two oscillations in the mean energy in Figs. Hand and is the typical duration 
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Figure 7. Turbulent lifetimes for different transitional Reynolds numbers, (a) Probability P(t) 
for a single trajectory to still be turbulent after a time t for six Reynolds numbers as indicated. 
Between 50 and 100 trajectories have been evaluated per Reynolds number. The distributions are 
well described by exponential distributions shown as straight lines, (b) Median r of the turbulent 
lifetimes as a function of Reynolds number. The median and the fluctuations increase rapidly 
with Reynolds number until the median reaches the cut-off lifetime of 2000 at Re — 2200. The 
error bars indicate the statistical uncertainty of the median. The straight dashed line shows the 
linear increase in lifetime expected due to purely non-normal linear dynamics. The inset shows 
the inverse median lifetime vs. Re and a linear fit, corresponding to a law r(Re) cx (Re c — Re) -1 , 
with Re c « 2250. 



of the regeneration cycle proposed in ^Hamilton et al\\l99& IWaleffell 19951 Il997t> . Thus, 
the Lyapunov exponent is small in absolute value, but fairly large on the intrinsic time 
scales of the turbulent dynamics. 

The extreme sensitivity of turbulent lifetimes to variations in initial conditions will 
make it next to impossible to prepare experimental perturbations sufficiently accurately 
to reproduce a run. A statistical analysis of, e.g., the distribution of lifetimes obtained by 
collecting data for several nearby initial conditions, should be more reliable. We therefore 
study P(t), the probability that a flow started with some initial condition will still be 
turbu lent after a time t. A related kind of statistics was proposed in ijDarbvshire fc Mullinl 
119951) : they analyzed the fraction p of initial conditions that remained turbulent over the 
time T it took for the perturbation to transit the distance between perturbation and 
detection. (Recall that we do not need to take this advection into account since our 
periodicity length is so short that the turbulence fills the entire volume.) In terms of our 
P(t), this probability is p = P(T). However, because of the unknown relation between 
the initial relaxation times in our simulations and the experiment we cannot compare 
values. 

For the lifetime statistics we used more than 50 initial amplitudes each for eight 
Reynolds numbers in the range from 1600 to 2200. The amplitudes were fairly high, 
w 1, in order to assure transition to turbulence at least for a short time interval. The 
results for different Reynolds numbers are shown in Fig. [7] 

The data support an exponential decay P(t) ~ exp(— et) for large times, with e the 
rate of escape from the turbulent state. An exponential distri bution of turbulent lifetimes 
is a characterist ic signature for escape from a chaotic saddle feadanoff fc TanJl984tlTel 
199 1t IOttlll993l). T hev have previously been seen in experiments (Bottin fc Chata ll998: 
Bott in et al\\l99£§ and nu merical studies on plan e Couette flow l|Schmiegellll999|) , as well 
as in Taylor-Couette flow l|Eckhardt et alm 002). They seem to be a generic feature of 
the transition in shear flows that are not dominated by linear instabilities. 
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Another measure that can be extracted from the lifetime distributions is the median, 
r = (ln2)/e. The median is the time up to which half the states have decayed. It is also 
well suited for numerical studies since it is not affected by the cut-off in integration time. 
The median increases rapidly with Reynolds number as shown in Fig. [TJd. Below Re = 
1800 the increase is mainly due to non-normal transient linear dynamics ijGrossmannl 
2000). When the Reynolds number is increased above Re — 2000 the median of the 
turbulent lifetimes as well as the fluctuations rise rapidly until the median reaches the 
cut-off lifetime of 2000 at Re ~ 2200. The inset in Fig.0shows the inverse median lifetime 
vs. Re and among the possible fits for these data the most satisfactory ones supports a 
divergence like (Re c — -Re) -1 , with Re c 2250. Such a divergence would be connected 
with a transition from the chaotic saddle for lower Re to an attractor for larger Re. 
However, it is k nown from other models that the lifetim es can increase rapidly without a 
true singularity llTellll99ltlCrutchfield Kanekolh 988^ . The question of whether we will 
arrive at a turbulent attractor cannot be answered here. But from the rapid increase it is 
clear that it will become an attractor for all practical purposes: for a setup with length 
1000 diameters the median will exceed this value for Reynolds number about 2130. For 
this observation time the fraction of states with shorter lifetimes drops from 80% of all 
initial conditions for Re = 2000 to 40% for Re = 2150 and to zero for Re=2200. These 
numbers are obtained from the numerical simulations and include the initial transient; 
numbers in other simulations and experiments may be different, but the rapid increase 
will remain the same. 



5. Conclusions 

An exponential distribution of lifetimes can be obtained for a constant probability of 
decaying. In view of the internal dynamics and the characteristic times connected with it, 
one can interpret this as the probability of decaying towards the laminar state at the end 
of each turbulent regeneration cycle. The conclusion from the exponential distribution of 
lifetimes is that this probability remains constant during the evolution, independent of 
the 'age' of the trajectory. 

It is worthwhile pointing out that the strong variation of the median lifetimes with Re 
is not reflected in the Lyapunov exponent, which increases only linearly with Re: this 
shows that the chaotic dynamics on the turbulent saddle and the escape from it are two 
different processes with independent characteristics. 

There are two predictions that should be accessible to experimental investigation: the 
behaviour near the boundary between laminar and turbulent and the exponential distri- 
bution of lifetimes. The increase in lifetime for increasing perturbation amplitude can be 
analyzed with a series of sensors along the pipe. The oscillations in Fig.E^ have a period 
of about 200 time units, which translates into 200 diameters. Thus, if a perturbation 
can be localized within a few radii, detectors that are further and further downstream 
might be able to reproduce some of the oscillations in Figure 4 as the perturbation ampli- 
tude is increased. Similarly, it might be possible to experimentally detect the Lyapunov 
exponents in the flow: with a Lyapunov exponent of about 0.07, a 5% uncertainty in 
the preparation of initial conditions will increase to 100% in a time of about 40 units, 
i.e. 40 diameters. Thus, placing detectors at spacings of several diameters apart should 
allow a determination of the separation between flow states starting from similar initial 
conditions. 

The second measurement is more easily performed: repeated runs with similar pertur- 
bations give the probability of finding a turbulent state that exists at least up to time t. 
The time can be varied again by placing detectors at different locations along the pipe. 
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As explained, the exponential distribution of lifetimes is a property of the turbulent state 
in the transition region and will not depend critically on the type of initial conditions 
used. Exact reproduction of an initial condition is hence not essential here since all initial 
conditions relax to the same turbulent state. The median lifetime increases rapidly with 
Reynolds number, from about 220 at Re = 1800 to 750 at Re — 2100. Given the ex- 
perimental limitations on the length of the pipe this region looks like a promising range 
for experiments. The mechanism studied here, the formation of a chaotic saddle, is fairly 
independent of the boundary conditions, but the quantitative characteristics may depend 
on it. Since the experiments will not have periodic boundary conditions in downstream 
direction and since the turbulence may be localized in puffs or slugs it will be interesting 
to see whether this will affect the lifetimes or its distribution. 

A detailed analysis of the distribution of lifetimes is a task which is still much better 
suited for experiments than for numerics: it is not necessary to prepare the initial distur- 
bance with highest precision, and the experimental observation time is of the order of 1 
minute for an experimental run, followed by about half an hour for the water to return 
to rest in the head-tank. On the other hand, a single trajectory in our computation takes 
up to 10 days. 

The picture that emerges from these data is that the transition to turbulence in pipe 
flow is connected with the formation of a chaotic saddle. The turbulent dynamics is 
chaotic, the lifetimes are exponentially distributed and the transition depends sensi- 
tively on the initial state. The same characteristics have been found in other shear flows 
without a linear instability as well llSchmiegel fc EckhardHll997t lEckhardt et aZJ 12002: 
iBottin et aZJIl99^: iFaisst fc Eckhardtll2000|) . The recent discovery of travelling waves in 
pipe flow for Reynolds numbers as low as 1250 l|Faisst fc Eckhardl^OOSUWedin fc Kerswelll 
2003) strengthens this picture by providing some of the states around which the network 
of homoclinic and heteroclinic connections that carry the chaotic saddle may form. 

Support by the German Science Foundation is greatfully acknowledged. 



REFERENCES 

Boberg, L. & Brosa, U. 1988 Onset of Turbulence in a Pipe. Z. Naturforsch. 43a, 697-726. 
Bottin, S. & Chate, H. 1998 Statistical analysis of the transition to turbulence in plane 

Couette flow. Eur. Phys. J. B 6, 143-155. 
Bottin, S., Daviaud, F., Manneville, P. & Dauchot, O. 1998 Discontinuous transition to 

spatiotemporal intermittency in plane Couette flow. Europhys. Lett. 43 (2), 171-176. 
Brosa, U. 1989 Turbulence without Strange Attractor. J. Stat. Phys. 55, 1303-1312. 
Canuto, C, Hussaini, M., Quarteroni, A. & Zang, T. 1988 Spectral Methods in Fluid 

Dynamics. Springer- Verlag. 
Crutchfield, J. & Kaneko, K. 1988 Are attractors relevant to turbulence? Phys. Rev. Lett. 

60, 2715-2718. 

Darbyshire, A. & Mullin, T. 1995 Transition to turbulence in constant-mass-flux pipe flow. 

J. Fluid Mech. 289, 83-114. 
Eckhardt, B., Faisst, H., Schmiegel, A. & Schumacher, J. 2002 Turbulence transition 

in shear flows. In Advances in turbulence IX (ed. I. Castro, P. Hancock & T. Thomas), p. 

701. Barcelona: CINME. 
Eckhardt, B. & Mersmann, A. 1999 Transition to turbulence in a shear flow. Phys. Rev. E 

60, 509-517. 

Eckhardt, B. & Yao, D. 1993 Local Lyapunov exponents. Physica D 65, 100-1008. 
Eggels, J., Unger, F., Weiss, M., Westerweel, J., Adrian, R., Friedrich, R. & Nieuw- 

STADT, F. 1994 Fully developed turbulent pipe flow: a comparison between direct numerical 

simulation and experiment. J. Fluid Mech. 268, 175-209. 



10 Holger Faisst and Bruno Eckhardt 

Faisst, H. & Eckhardt, B. 2000 Transition from the Couette- Taylor system to the plane 

Couette system. Phys. Rev. E 61, 7227. 
Faisst, H. & Eckhardt, B. 2003 Traveling waves in pipe flow. Phys. Rev. Lett. 91, 224502. 
GROSSMANN, S. 2000 The onset of shear flow turbulence. Rev. Mod. Phys. 72, 603-618. 
Hamilton, J., Kim, J. & Waleffe, F. 1995 Regeneration mechanisms of near-wall turbulence 

structures. J. Fluid Mech. 287, 317-348. 
Hof, B., Juel, A. & Mullin, T. 2003 Scaling of the turbulence transition threshold in a pipe. 

Phys. Rev. Lett. 91, 244502. 
Kadanoff, L. & Tang, C. 1984 Escape from strange repellers. Proc. Natl. Acad. Sci. USA 81, 

1276. 

Ma, B., van Doorne, C, Zhang, Z. & Nieuwstadt, F. 1999 On the spatial evolution of a 
wall-imposed periodic disturbance in pipe Poiseuille flow at Re — 3000. Part 1. Subcritical 
disturbance. J. Fluid Mech. 398, 181-224. 

Meseguer, A. & Trefethen, L. 2003 Linearized pipe flow to Reynolds number 10.000.000. 
J. Comput. Phys. 186, 178-197. 

Ott, E. 1993 Chaos in Dynamical Systems. Cambridge University Press. 

Priymak, V. & Miyazaki, T. 1998 Accurate Navier-Stokes investigation of transitional and 

turbulent flows in a circular pipe. J. Comp. Phys. 142, 370-411. 
Quadrio, M. & Sibilla, S. 2000 Numerical simulation of turbulent flow in a pipe oscillating 

around its axis. J. Fluid Mech. 424, 217-241. 
Reynolds, O. 1883 An experimental investigation of the circumstances which determine 

whether the motion of water shall be direct or sinuous and the law of resistance in parallel 

channels. Phil. Trans. R. Soc. 174, 935-982. 
Rubin, Y., Wygnanski, I. & Haritonidis, J. 1980 Further observations on transition in a 

pipe. In Laminar- Turbulent Transition (ed. R. Eppler & F. Hussein), pp. 19-26. Springer. 
Schmid, P. & Henningson, D. 1994 Optimal energy density growth in Hagen-Poiseuille flow. 

J. Fluid Mech. 277, 197-225. 
Schmiegel, A. 1999 Transition to turbulence in linearly stable shear flows. PhD thesis, Philipps- 

Universitat Marburg. 

Schmiegel, A. & Eckhardt, B. 1997 Fractal stability border in plane Couette flow. Phys. 

Rev. Lett. 79 (26), 5250-5253. 
Schmiegel, A. & Eckhardt, B. 2000 Persistent turbulence in annealed plane Couette flow. 

Europhys. Lett. 51, 395-400. 
Tel, T. 1991 Transient chaos. In Directions in Chaos (ed. H. Bai-Lin, D. Feng & J. Yuan), , 

vol. 3. World Scientific, Singapore. 
Trefethen, L., Chapman, S., Henningson, D., Meseguer, A., Mullin, T. & Nieuw- 
stadt, F. 2000 Threshold amplitudes for transition to turbulence in a pipe, http://xxx.uni- 

augsburg.de/abs/physics/0007092 . 
Trefethen, L., Trefethen, A., Reddy, S. & Driscol, T. 1993 Hydrodynamics stability 

without eigenvalues. Science 261, 578-584. 
Waleffe, F. 1995 Transition in shear flows. Nonlinear normality versus non-normal linearity. 

Phys. Fluids 7, 3060-3066. 
Waleffe, F. 1997 On a self-sustaining process in shear flows. Phys. Fluids 9, 883-900. 
Wedin, H. & Kerswell, R. 2003 Exact coherent structures in pipe flow: travelling wave 

solutions, submitted . 

Wygnanski, I. & Champagne, F. 1973 On transition in a pipe. Part 1. The origin of puffs 
and slugs and the flow in a turbulent slug. J. Fluid Mech. 59, 281-335. 

Wygnanski, I., Sokolov, M. & Friedman, D. 1975 On transition in a pipe. Part 2. The 
equilibrium puff. J. Fluid Mech. 69, 283-304. 

Zikanov, O. 1996 On the instability of pipe Poiseuille flow. Phys. Fluids 8 (11), 2923-2932. 



